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The electromagnetic properties of the S'L r (3)-fiavor baryon decuplet are exam- 
ined within a lattice simulation of quenched QCD. Electric charge radii, magnetic 
moments, and magnetic radii are extracted from the EO and Ml form factors. Pre- 
liminary results for the E2 and M3 moments are presented giving the first model 
independent insight to the shape of the quark distribution in the baryon ground state. 
As in our octet baryon analysis, the lattice results give evidence of spin-dependent 
forces and mass effects in the electromagnetic properties. The quark charge distribu- 
tion radii indicate these effects act in opposing directions. Some baryon dependence 
of the effective quark magnetic moments is seen. However, this dependence in decu- 
plet baryons is more subtle than that for octet baryons. Of particular interest are the 
lattice predictions for the magnetic moments of and A ++ for which new recent 
experimental measurements are available. The lattice prediction of the A ++ /p ratio 
appears larger than the experimental ratio, while the lattice prediction for the 0~ /p 
magnetic moment ratio is in good agreement with the experimental ratio. 
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I. INTRODUCTION 



Knowledge of the quark substructure of baryons is largely based on experiment and 
model dependent descriptions of quark-gluon interactions. While most models are "QCD 
inspired", they all suffer from an approximate treatment of the complex nonperturbative 
long distance interactions. Often these interactions are approximated by a simple potential. 
Quark spin dependent interactions are approximated by a single gluon or effective pion ex- 
change. In many models the subtleties of current quark and gluon interactions are lumped 
into an effective or constituent quark mass. Skyrme models offer an alternative description 
of hadronic phenomena. However, the foundation of the model is in the 1/N C expansion 
of QCD. To learn the true nature of nonperturbative QCD and, we hope, hadronic phe- 
nomena, it is necessary to calculate directly with the QCD Lagrangian in a manner which 
fully accounts for nonperturbative interactions. The most successful, reliable and promising 
approach currently available is that of numerical simulations of QCD. 

Through studies of hadronic electromagnetic form factors, the lattice gauge approach to 
QCD has proved to be a valuable tool in revealing the underlying quark substructure of 
hadrons ^M- Early calculations focused on the pion electric form factor with SU{2) color 
0-^1] and later with SU (3) color Calculations of the proton electric form factor followed 

0. Electromagnetic form factors of n, p and N were calculated from which magnetic 
moments and electric charge radii were extracted. Our analysis of the entire baryon flavor- 
octet followed |TJ in which electromagnetic properties were reported for both baryons and 
the quark sector contributions. This examination of octet baryons exposed the presence of 
spin-dependent forces and center-of-mass effects in the underlying quark dynamics. These 
effects give rise to large variations in the quark contributions to baryon magnetic moments 
which were not anticipated by model calculations. Recently the q 2 dependence of the nucleon 
electromagnetic form factors was examined using a method which characterizes one of the 
nucleon interpolating fields as a zero momentum secondary source [|10| . 

Other studies of hadron structure in lattice QCD have been pursued through an exami- 
nation of current overlap distribution functions [] I -13|1 and Bethe-Salpeter amplitudes [13 



Form factors have a number of advantages over Bethe-Salpeter amplitudes, in that they are 
gauge invariant, path independent and allow the extraction of quark distribution radii rela- 
tive to the system center of mass. In contrast, Bethe-Salpeter amplitudes only give relative 
quark separations and subtle dynamical effects can remain hidden. 

In this paper we continue our study of hadron structure and present the first lattice 
QCD calculation of the electromagnetic form factors of 5'?7(3)-flavor-decuplet baryons. An 
analysis of electromagnetic transition moments will follow in a subsequent paper fl4j| . The 
technique for extracting the four form factors associated with spin-3/2 baryons from the 
electromagnetic current matrix elements has been outlined in Ref. [I5||. From these form 



factors we will determine magnetic moments, electric radii and magnetic radii and present 
preliminary results for the higher-order E2 and M3 moments. 

On the lattice, decuplet baryons are stable as a result of the unphysically large quark 
masses that are used in present calculations. Decay to a pion and an octet baryon is 
forbidden by energy conservation. However, stability of decuplet baryons is common to most 
hadronic models. In this sense, lattice results provide a new forum in which the strengths 
and weaknesses of various models may be identified. The lattice results also provide access to 
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observables not readily available with present experiments such as the higher-order multipole 
moments of the Q~ which is stable to strong interactions. 

An examination of decuplet baryon structure in lattice QCD enables one to study new 
aspects of nonperturbative quark-gluon dynamics. In analyzing the results we make com- 
parisons within the baryon decuplet and with the octet results which provide insight into the 
spin dependence of quark interactions. The E2 and M3 moments accessible in spin-3/2 sys- 
tems provide a preliminary glimpse at the shape of the decuplet baryon ground state. These 
higher-order moments also have the potential to discriminate between model dependent de- 
scriptions of hadronic phenomena. For example, a vanishing E2 moment would cast serious 
doubt on hedgehog Skyrmion descriptions of baryons. To put our results into perspective, 
we compare our calculations with experimental measurements where available, with recent 
quark, bag and Skyrme model calculations, and with QCD sum rule calculations. 

The format of our paper is as follows. Interpolating fields are discussed in section II 
A. Correlation functions at the quark level are discussed in section II B. Two-point and 
three-point correlation functions at the hadronic level are discussed in sections II C and II D 
respectively. Lattice techniques are discussed briefly in section II E. Decuplet baryon masses 
are reported in section III A. Correlation function ratios used in extracting the multipole 
form factors are illustrated in section III B. Our findings for the electromagnetic properties 
of the four multipole form factors are presented in sections III C through III F. Finally, 
section IV provides an overview of our results and a discussion of future investigations. 



II. THEORETICAL FORMALISM 



A. Interpolating fields 



The commonly used interpolating field for exciting the A ++ resonance from the QCD 
vacuum takes the long established ||T6| , |T7f form of 



X^ {x) 



abc f„,Ta 



u Ia (x)C^u b (x))u c (x). 



(2.1) 



Unless otherwise noted, we follow the notation of Sakurai The Dirac gamma matrices 
are Hermitian and satisfy {7^,7^} = 2 5^, with a^ u = — [7^,7^]. G = 7472 is the charge 
conjugation matrix, a, b, c are color indices, u(x) is a w-quark field, and the superscript T 
denotes transpose. The generalization of this interpolating field for the A + composed of two 
u quarks and one d quark has the form 



xt + ( x ) 



1 
73 



ubc 



2 (u Ta (x)C lpL d\x)) u c (x) 
+ (u Ta (x)C^u b (xj) d c {x) 



(2.2) 



Other decuplet baryon interpolating fields are obtained with the appropriate substitutions of 
u(x), d(x) — > u(x), d(x) or s(x). The interpolating field for S*° is given by the symmetric 
generalization 



- e 



abc 



u Ta (x)C^d b (x)) s c {x) 



.3 



+ (d Ta (x)C-y^s b (x)) u c {x) 
+ (s Ta (x)C 7 X(^)) d c {x) 
The 5'C/(2)-isospin symmetry relationship for S* form factors 



i*0 



(2.3) 



(2.4) 



may be easily seen in the S*° interpolating field by noting 

e abc (sT*( x )C lllU \x)) d C {x) = 



^abc 



(u Ta (x)C lfM s\x)) d c (x) 



(2.5) 



B. Correlation functions 

Two-point correlation functions at the quark level are obtained through the standard 
procedure of contracting out pairs of quark fields. Consider the A + correlation function at 
the quark level. 



T(xf(x),xf (0) 



1 

3 

+ 
+ 
+ 



e abc e 



4S° 

V 

AS" 



45; 



2S7 
25" 



2S a d a 



lv cs T d bb 'c ltl s c : 
lv cs T u bb> c 1 ,s c : 



(2.6) 



tr 
tr 
tr 



lu CS T d bb 'C 1 ,S c u 
lu CS T u bb 'C 1 ,S^']} 



where the quark-propag ator S™ = T (u a (x), u a ' (0)j and similarly for other quark flavors. 
5'L r (3)-flavor symmetry is clearly displayed in this equation. 

The corresponding connected three-point function may be constructed by replacing each 
of the three propagators S, one at a time, by S denoting the propagation of a quark in 
the presence of the electromagnetic current. The interesting result obtained from (|2.6| ) is 
that the electromagnetic form factors of the neutral charge decuplet baryons vanish under 
S'?7(3)-flavor symmetry. Introduction of the more massive strange quark, as in S*° or S*°, 
violates 5'[/(3)-flavor symmetry and allows a nontrivial result. However, under SU(2)- 
isospin symmetry the electromagnetic form factors of A vanish. A non-zero value for the 
magnetic moment of the A resonance reflects differences in the u- and <i-quark masses and 
contributions from the quark sea through disconnected quark loops interacting with the 
electromagnetic current. 
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C. Two-point Green functions 



In this and the following subsection discussing correlation functions at the hadronic level, 
the Dirac representation of the 7-matrices as defined in Itzykson and Zuber |19| is used to 
facilitate calculations of the 7-matrix algebra. It is then a simple task to account for the 
differences in 7-matrix and metric definitions in reporting the final results using Sakurai's 
notation. 

The extraction of baryon mass and electromagnetic form factors proceeds through the 
calculation of the ensemble average (denoted (•■■)) of two and three-point Green functions. 
The two-point function is defined as 

(0;p;r))= 

£ e-^r^ ( Q I T (x a a (x)xm) I ) • (2-7) 

X 

Here Q represents the QCD vacuum, T is a 4 x 4 matrix in Dirac space and a, (3 are 
Dirac indices. The subscripts cr, r are the Lorentz indices of the spin-3/2 interpolat- 
ing fields. At the hadronic level one proceeds by inserting a complete set of states 
I B,p, s ) and defining 

( n I X a{0) \B,p,s)= Z B J ^- u a (p, s), (2.8) 

V 

where Zb represents the coupling strength of x(0) to baryon B. Momentum is denoted by 
p, spin by s, and u a (p,s) is a spin-vector in the Rarita-Schwinger formalism [p0[| . E p = 
^Jp 1 + M 2 and Dirac indices have been suppressed. Using the Rarita-Schwinger spin sum, 

P + M \ 1 2p a p T Pa^r-l 

2M { 9aT 3 7a7r 3M 2 + 3M 

our usual definitions for T, 

°) r =-( I 

2 1 ' 4 2 1 



L MP, s)u T (p, s) = — — g aT - - 7 ,7r - ^rry + ^rj , (2-9) 



ri = ; °n J „ ; r 4 = ^( J n ^ , (2.10) 



and p = (p, 0, 0), the large Euclidean time limit of the two point function takes the form 

( G B J{t-p, r 4 ) ) = zl^e-^tr [ r 4 K T ] , (2.11) 



where 



2 |p1 2 [E p + M B \„ t 



Gtfmr A ) ) =Zy-^\ ^— — e-^, (2.12a) 



3 Ml V 2E P 



^(t;^r 4 ))^^g(^je--, (2.12b) 



( Gg B (t;p,T,) ) =Z\\ {^^j e-*», (2.12c) 
( Gg B (t;p,T 4 ) > = Z\\ (^^) e-^- (2.12d) 
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In determining the appropriate forms suitable for calculations using Sakurai's conventions 
the definitions of the 7-matrices used in the interpolating fields are taken into account. Since 
the nonvanishing terms of are diagonal in a and r, the 7-matrices are paired with their 
Hermitian conjugates. Since the 7-matrix notations differ only by factors of i and —1, there 
are no required alterations for calculations using the notation of Sakurai. 

Equations (|2.12a|) through fl2.12d| ) provide four correlation functions from which a baryon 
mass may be extracted. All baryon masses extracted from the different selections of Lorentz 
indices agree within statistical uncertainties. The combination providing the smallest sta- 
tistical fluctuations is ( G^it; p, T 4 ) + Gf^(t; p, T 4 ) ) and these results are presented in 
section III. 

It should be noted that the spin-3/2 interpolating field also has overlap with spin-1/2 
baryons. For the A baryons and Q~ this poses no problem as these baryons are the lowest 
lying baryons in the mass spectrum having the appropriate isospin and strangeness quan- 
tum numbers. However, £* and S* correlation functions may have lower lying spin-1/2 
components and therefore it is desirable to use the spin-3/2 projection operator [^1[ 



p % 2 (p) = 9n» ~ 3 liflv ~ ^2" ■ P In Vv + Pn J v 7 ■ p) • (2-13) 

However, to use this operator, one must calculate the full 4 4 matrix in Dirac and Lorentz 
spaces of G^(t;p,T) which exceeds our current analysis of 4 Lorentz terms and 2 Dirac 
terms by a factor of 32. 

QCD sum rule investigations of S* and S* hyperons suggest that the spin-1/2 component 
of the spin-3/2 interpolating field is small relative to the spin-3/2 component |22|| . How- 
ever, the analysis does not determine whether the spin-1/2 component lies above or below 
the lowest lying spin-3/2 state. Our lattice results for baryon two-point functions give no 
indication of a low- lying spin-1/2 component being excited by the spin-3/2 interpolating 
fields. 



D. Three-point functions and multipole form factors 



Here we begin with a brief overview of the results of Ref. ||15|| , where the multipole form 
factors are defined in terms of the covariant vertex functions and in terms of the current 
matrix elements. The Dirac representation of the 7-matrices as defined in Itzykson and 
Zuber continues to be used to facilitate calculations of the 7-matrix algebra. Finally, the 
results are reported in Sakurai's notation in a form suitable for calculation in lattice field 
theory. 

The electromagnetic current matrix element for spin-3/2 particles may be written as 



( p',s' I /(O) \p,s)~- 



M 2 

" B u a (p',s')O a ^up(p,s). 



\ EpEj 



(2.14) 



Here p, p' denote momenta, s s' spins, and u a (p, s) is a Rarita-Schwinger spin- vector. The 
following Lorentz covariant form for the tensor 



°^ - -<* + ^} " A + 2^} 



(2.15) 
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where P = p' + p, q = p' — p and Mb is the mass of the baryon, satisfies the standard 
requirements of invariance under time reversal, parity, G-parity and gauge transformations. 
The parameters a 1; a 2 , c\ and c 2 are independent covariant vertex function coefficients which 
are related to the multipole form factors. 

The multipole expansion of the electromagnetic current matrix element, defined in terms 
of angular momentum recoupling algebra, has the following form 



(p',s'\f(0)\p,s) = 

A (l s '\ G Eo(q 2 ) + 2V5rG E2 (q 2 ) [E< a > x|x q]^ 



(o) , 3 



(2.16a) 



(p',s' | J(0) \p,s )= 

( \s' | \G EQ {q 2 ) + 2V5rG E2 (q 2 ) [e™ x [g x gp) 

+ *{^Wg 2 )£« + 3rG A/3 (g 2 ) [E<» x [qxqf^ 



(0) 



P 



x Q I 2 s ) 



(2.16b) 



where r = — q 2 /{2Mb) 2 (> 0), and P and q are unit vectors. A = y/1 + r in the laboratory 
frame (p = 0) and A = 1 in the baryon Breit frame (P = p' + p = 0). The spin matrix 
elements are defined by Clebsch-Gordan coefficients, 



(^'|SW|^0=v / 15(^lm 



- f-s'2m 
6 V2 



3 3 
2 X 2 S 

his 
2 2 



--\ - (~s'3m 



6 V 3 V2 



3 q 3 
-3-s 
2 2 



(2.17a) 
(2.17b) 

(2.17c) 
(2.17d) 



where the Condon and Shortley phase convention has been used [23 . 

The multipole form factors are defined in terms of the covariant vertex function coeffi- 



cients ai, a 2 , C\ and c 2 through the following Lorentz invariant expressions |L5 



;i + \r) {ax + (1 + r)a 2 } - ^r(l + r) { Cl + (1 + r)c 2 } 



£mi(<? 2 ) 



{a x + (1 + r)a 2 } - ^(1 + r) {ci + (1 + r)c 2 } , 

4 2 
(1 + gT")ai - + t-)ci , 

oi - -(1 + r)ci . 



(2.18a) 
(2.18b) 
(2.18c) 
(2.18d) 



The multipole form factors Geo, Ge2, Qmi and Qm3 are referred to as charge (E0), electric- 
quadrupole (E2), magnetic-dipole (Ml) and magnetic-octupole (M3) multipole form fac- 
tors, respectively. 
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In a manner similar to that for the two-point function, the three-point Green function 
for the electromagnetic current is defined as 

(G^ B (t 2 ,t i; p',p;T)) = 

j2 e -iP-*\+iv'-?>*re<* < n | t fax^ffayg (o)) | n > . (2.19) 

Once again, the subscripts a, r are the Lorentz indices of the spin-3/2 interpolating fields. 
For large Euclidean time separations t 2 — t\ » 1 and t\ » 1 the three-point function at 
the hadronic level takes the limit 

( G^ B (t2,h;p',p;T) >= £ e - E P '(^) e -^ 

s,s' 

r? a ( n | x« I p',s' )(p',s' \r\p,s)(p,s\xP\tt). (2.20) 

where the matrix element of the electromagnetic current is defined in ( |2.14| ), and the matrix 
elements of the interpolating fields are defined by ( |2.8| ). 

In Ref. |]T3| it was noted that the time dependence of the three-point function may be 
eliminated by use of the two-point functions. However the appropriate combination of the 
two-point function Lorentz indices was not specified. Maintaining the lattice Ward identity, 
which guarantees the lattice electric form factor reproduces the total charge of the baryon 
at q 2 = 0, provides an indispensable guide to the optimum ratio of Green functions. The 
preferred ratio of two- and three-point Green functions is 



R^At 2 ,t i; p',p;T) = 

G*r B (t2,h;p',p;r) }(G B r B (t2,t 1 ;-p,-p';T) }\ 1/2 ^ 



(G B T B (t 2 ;p';T 4 ))(G^(t 2 ;-p;T A )} 

' E p + M \ 1/2 f E p , + M \ 1/2 - 
2E P J \ 2E p/ J 



R a \(p',p;T), (2.22) 



where we have defined the reduced ratio R a fl T (p',p; T). Note that there is no implied sum 
over cr and r in (|2.21|). 



Using our standard definitions for T given in ( |2.10|) and the Rarita-Schwinger spin sum 
of (|2.9| ), the multipole form factors may be isolated and extracted. The appropriate com- 
binations of R a ^ T {p',p\ T) suitable for calculations employing the 7-matrix and metric con- 
ventions of Sakurai are 

GEo(q 2 ) = I (Ri\(qi, 0; T 4 ) + R 2 \(q u 0; T 4 ) + R 3 \(qi, 0; T 4 )) , (2.23a) 
Q E2 {q 2 ) = 2 M(E + M) (i?i 4 i(<fi, 0; T 4 ) + R 2 \(q u 0; T 4 ) - 2R 3 \(q u 0; T 4 )) , 

ml 

(2.23b) 

QmM 2 ) = - f^i 3 i(gi, 0; r 2 ) + i? 2 3 2 (gi, 0; r 2 ) + R 3 3 3 (q u o ; r 2 )) , (2.23c) 

5 |gi| v > 
Quz{q 2 ) = _ 4 M( ^+ M)2 fSi 3 i(gi,0;r 2 ) + E 2 3 2 (gi, 0; T 2 ) - | i? 3 3 3 (gi, 0; T 2 ^ 

(2.23d) 



qi\ s V 2 



s 



where 

/ 3 „/ I 



qi 



,(2) 



(g, 0,0). Equation ( |2.23b| ) for ^2 isolates the spin matrix element 



)s ). Smaller statistical uncertainties may be obtained by using the sym- 
metry -R2 4 2(<fi, 0; T 4 ) = i? 3 4 3 (g!, 0; T 4 ). This symmetry is used in the following analysis to 
eliminate the term -R2 4 2(<?i, 0; T 4 ) in the E2 form factor. The signs of the El form factors 
remain unchanged. 

It is worth noting at this point that the form of the tensor O a ^ in ( 2.15 ) is not unique. 
While equations (|2.23a|) through (|2.23d|) are derived using fl2.15|) these results are more 
generally applicable. For example, the authors of Ref. |24j] employ the form 



q a g (3 
2M| 



{ F 3 >f + 



-i(T ta, q v 



2M, 



B 



(2.24) 



Provided one includes a factor of v^6 to account for their normalization of the M3 form 
factor, the results of (|2.23aj) through fl2.23d| ) isolate the multipole form factors. 



E. Lattice techniques 

Here we briefly summarize the lattice techniques used in the following calculations. Ad- 
ditional details may be found in Ref. [[J . 

We use Wilson's formulation for both the gauge and fermionic action. Si7(2)-isospin 
symmetry is enforced by equating the Wilson hopping parameters k u = kj, = k. We select 
three values of k, which we denote K\ = 0.152, n 2 = 0.154 and K3 = 0.156, and extrapolate 
the u-d quark sector to the chiral limit. To account for the relatively heavy strange quark 
we fix k s = Ki, the smallest of the three values of k considered. This allows an accept- 
able extrapolation of the light quarks to the chiral limit through values of quark mass less 
than or equal to the strange quark mass. Our calculations of octet and decuplet baryon 
masses indicate that this selection for k s gives a reasonable description of the strange quark 
dynamics. 

The conserved electromagnetic current is derived from the fermionic action by the 
Noether procedure. The lattice Ward identity guarantees the lattice electric form factor 
reproduces the total baryon charge at q 2 = 0. The quark propagators coupled with fixed 
momentum q\ = (q, 0, 0) to are calculated using the sequential source technique (SST) 

To minimize noise in the Green functions, we exploit the parity symmetry of the correla- 
tion functions, and the equal weighting of {U} and {U*} gauge configurations in the lattice 
action. Defining sp as 

G(p',p,q;T) = Sp G(-p>,-p,-q;T), s P = ±1, (2.25) 
and sc as 

T = s c (CTC- 1 )\ s c = ±l, (2.26) 
where C = C75, the correlation functions are real provided 
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s c = s P . (2.27) 

This condition is satisfied with the selections for T indicated in ( |2.23a| ), though d2.23d| ). 
While this approach requires an extra matrix inversion to determine an additional SST 
propagator with momentum —q\, inclusion of both {U} and {U*} configurations in the 
calculation of the correlation functions provides an unbiased estimate of the ensemble average 



properties which has substantially reduced fluctuations p8fl . 

Twenty-eight quenched gauge configurations are generated by the Cabibbo-Marinari [23 
pseudo-heat-bath method on a 24 x 12 x 12 x 24 periodic lattice at (3 — 5.9. Dirichlet 
boundary conditions are used for fermions in the time direction. Configurations are selected 



after 5000 thermalization sweeps from a cold start, and every 1000 sweeps thereafter |30 
Time slices are labeled from 1 to 24, with the 5-function source at t = 4. A symmetric 
combination of the current (j M (xi — p) + j fJ- (xi))/2 is centered at time slice t\ = 12. The 
following calculations are done in the lab frame p = 0, p' = q\ = (27r/24, 0, 0), the minimum 
nonzero momentum available on our lattice. The spatial direction of the electromagnetic 
current is chosen in the z-direction. As discussed in section III B and as in Ref. [|l|], the 
reported masses and form factors are determined by fitting the correlation functions in time 
slices 16 through 20 inclusive. 

Statistical uncertainties are calculated in a third-order, single elimination jackknife 



31.32]. A third order jackknife provides uncertainty estimates for the correlation functions, 



fits to the correlation functions, and quantities extrapolated to the chiral limit. 

III. RESULTS 
A. Baryon mass 

Baryon masses are determined by fitting the Euclidean time evolution of baryon two- 
point functions at three values of k with the strange quark fixed at the smallest value of 
k. The light u and d quarks are extrapolated to n cr where an extrapolation of the squared 
pion mass vanishes. The nucleon mass is used to define the lattice spacing a = 0.128(11) 
fm, a -1 = 1.54(13) GeV. Figure [I] displays the masses of the decuplet baryons obtained 
from the combination of two-point functions + Gf^ . Octet baryon masses are also 



displayed for reference. Experimental masses |53[] are indicated by horizontal lines. With 
the exception of the N-A splitting, the baryon masses are reasonably reproduced. Table 
| summarizes the lattice baryon masses at the three values of k considered as well as the 
extrapolated masses. The momentum transfer is relatively insensitive to the baryon mass. 
For all baryons, q 2 a 2 = 0.067 ± 0.001 with the larger values corresponding to the heavier 
baryons. 

B. Correlation function ratios 



Here we present a few of the correlation function ratios of Q2.21 ) which are combined to 



isolate the multipole form factors. Figure |2| displays the quark contributions to the electric 
charge form factor of £* hyperons at k u = k 3 , the lightest quark masses considered. The 
ci-quark contribution is related to the u-quark result by a ratio of the quark charges. Electric 
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charge form factors for the three charge states of the S* hyperons are obtained through the 
appropriate sums of quark contributions. Similar results are seen for other decuplet baryons. 

The form factors are determined by fitting the correlation function ratio sum by a hor- 
izontal line for times t 2 » 12 and t 2 « 24. Fortunately there is a rather broad plateau 
region where the electric form factor may be reliably determined. We consider fits of the 
correlation functions from time slice 15 through 21 in intervals including 4 to 7 points. The 
results are selected from these 10 fits based on the flatness of the correlation functions and 
the statistical uncertainties. As in the octet baryon analysis, it is found that fits of the 
5 points in the time slice interval 16 to 20 provide the optimum balance between these 
systematic and statistical uncertainties. 

Figure ||| shows a similar plot for the quark contributions to the magnetic-dipole form 
factor of S* hyperons. Once again k u = k 3 with n s fixed at K±. 

Figures (| and |5] display the E2 and M3 ratio sums for A baryons in comparison with E0 
and Ml ratio sums respectively. Only the ratio terms appearing inside the parentheses of 
( p.23a| ) through (|2.23d| ) are included here. The higher-order moments are dependent upon 
subtle differences in the ratio contributions, R a ^ T . Fortunately the statistical uncertainties of 
individual correlation function ratios are correlated and some of the uncertainty is canceled 
in the higher-order moment sums. However, statistical uncertainties in the E2 and M3 
moments remain large. 

The El result is small relative to the E0 ratio sum. In the interest of placing an upper 
limit on the magnitude of the E2 form factor we will consider fits of the E2 correlation 
functions in the standard interval of t 2 = 16 through 20. The M3 moment is finite and 
remains finite for all three values of k considered. This is not entirely surprising as pure 
Ml multipole moments are generally not seen in extended physical systems. Unfortunately, 
a clear plateau region is not seen. To determine the scale of the M3 moment we will fit 
the correlation function ratio sum in the standard interval of t 2 = 16 through 20 with the 
understanding that the result may be an underestimate of the M3 moment. 



C. Electric charge form factors 

Electric charge form factors are extracted from fits of the ratios in ( |2.23a| ). Table [II 
reports the form factors for the decuplet baryons at the three values of k considered along 
with the extrapolated values. Table |T| gives the form factors for individual quarks of unit 
charge within decuplet baryons. Note that isospin symmetry equates ci-quark properties to 
u-quark properties when normalized to unit charge. 

We use the standard small q 2 expansion of the Fourier transform of a spherical charge 
distribution to extract the electric charge radius. 

( r 2 > = -^ 2 g E (q 2 ) \ q2=0 . (3.1) 

We have two points describing the function QeIq 2 )- Namely (?e(0), the total charge of the 
baryon, and Qe{q 2 ) evaluated at the smallest finite value of q 2 available on our lattice. To 
evaluate the derivative of (|3.1|) we use a dipole form for the function Q e (q 2 )- The dipole 
result for the radius is 



11 



(r 2 ) 
Ge(0) 




(3.2) 



To assess the sensitivity of our results on the dipole approximation we also consider a 
monopole form. 



In the figures and tables of this section we quote the quantity J ( r 2 ) /Qe(Q) which gives 
the radius of baryons and quark distributions with unit charge. In all cases the sign of ( r 2 ) 
is the same as the charge of the baryon or quark. The difference of the radii extracted in 
the dipole and monopole approximations is small relative to the statistical uncertainties in 
the radii. We refer to the dipole results in the following discussion and figures. 

The charge radii of decuplet baryons are not determined by experiment so our goal here is 
to compare decuplet and octet results to see if there is some indication of how spin- dependent 
interactions affect the electromagnetic structure. Of course the spin-dependent interactions 
depend on the quark masses and our calculations are done at masses which are not yet 
small. However the extrapolation of charge radii into the physical (light) quark region is 
problematic since charge radii become infinite in the exact chiral limit ||34|| . To reproduce 



this behavior in the present calculation would require some model dependent theoretical 
input which is not available. 

To estimate the charge radii in the physical mass region and allow a comparison with 
model calculations which often neglect couplings to the continuum, we extrapolate in l/« (or 
equivalently in m 2 ) to n cr . Although this prescription involves some systematic uncertainty 
p5| , it gives a better picture of what is happening at physical quark masses than, for example, 
using only results calculated at our largest k value. However it should be noted that the 
qualitative statements and conclusions are the same whether we extrapolate or not. 

Figure |6] displays the lattice results for the electric charge radii normalized to the proton 
charge radius for the charged members of the baryon octet and decuplet. We find a similar 
pattern of quark mass effects in the octet and decuplet with a A + charge radius essentially 
equal to that of the proton. Note that the charge radii of A ++ and A~ are equal to that 



of A + when the radii are normalized by the total charge. Table [TV] summarizes the lattice 
calculations. 

Some insight into quark dynamics may be obtained by examining the behavior of the 
quark distribution radii in hyperons as the u and d quarks become lighter. Consider for 
example, the distribution radii of u and s quarks within S + as the u quarks are extrapolated 
to n cr . Figures [F] and [8] display the quark distribution radii in lattice units (LU) for the octet 
S + and the decuplet S* + respectively. 

In both cases, the radius of the w-quark distribution increases as the w-quark becomes 
lighter (smaller k _1 ), as expected. However there is a difference in the behavior of the s- 
quark distributions in the octet E + and the decuplet In the case of octet E + , there is a 
small, but significant, trend for the s-quark distribution radius to decrease as the n-quarks 
become lighter ||36|| . This effect was understood as a shifting of the center of mass of the u-s 
system towards the strange quark which becomes relatively heavier as the w-quarks become 
lighter. In contrast, the s-quark distribution radius in the decuplet S* + has little dependence 
on the M-quark mass. 

This difference between octet and decuplet behavior has a natural explanation in the 
hyperfine interaction term of the one-gluon-exchange potential. The hyperfine force is re- 
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pulsive for quarks interacting in spin-triplet states. The strength of the interaction increases 
with decreasing -u-quark mass and provides a mechanism to counteract the center-of-mass 
effect. 

Lattice results for the neutron electric form factor confirm a repulsive force between 
quarks with their spins aligned. Furthermore, the magnetic form factors indicate that in 
the octet S + (decuplet the singly represented s quark has, on average, its spin anti- 

aligned (aligned) with that of the doubly represented u quarks. Hence the spin alignments of 
u and s quarks seen in the lattice results are in qualitative agreement with those anticipated 
from S , f/(6)-spin-flavor symmetry. Note however, the lattice dynamics do not impose this 
symmetry. The quark contributions to octet baryon magnetic moments differ significantly 
from the SU(6) predictions. 

Figure |9] shows charge distribution radii for quarks of unit charge in decuplet baryons. 
The analogous graph for octet baryons is indicated in figure [10]. It should be noted that 
although the uncertainty regions of the radii for different quarks overlap it does not neces- 
sarily mean that the u quark distribution radius in n, for example, may be larger than the 
u quark radius in S°. The uncertainties are highly correlated between the two results and a 
calculation of the difference of the radii indicates u~ is larger by 0.50 +J J LU. Other octet 
baryon quark pairs such as u^ — u p differ from zero by at least one standard deviation. In the 
octet, therefore, some nontrivial baryon dependence of the quark distributions does occur. 
In contrast, significant baryon dependence of the quark distributions in decuplet baryons is 
not observed. The lattice results are summarized in table [V]. 

Figure |TT] displays ratios of octet and decuplet quark charge distribution radii. With the 
exception of the singly represented octet quarks u n and ss which showed some baryon de- 
pendence in their charge distributions, the remaining quark distribution radii are unaffected 
by differences between octet and decuplet spin-flavor symmetry. 

Figure |Tl| also gives an understanding as to why the proton and A + have similar charge 
radii. Naively one might expect A + to be larger due to additional spin-dependent repulsion 
between the quarks. However, the dominant effect of the additional spin repulsion is to 
enhance the charge distribution radius of the d quark in A + relative to that in p. This 
counteracts any possible enhancement of u-quark radii in A + . 

Some attention has been given to model calculations of the Q~ charge radius. The model 
predictions vary over a wide range. Figure [12] illustrates results for the Q~ jp charge radius 
ratio. The results include the lattice (Latt.) results of this analysis, a calculation based 
on the relativised quark model (Q.M.) of Ref. [0, MIT bag (MIT) and cloudy bag (C.B.) 
models |38| and a bound state approach Skyrme (Skyr.) model calculation [j39| . Only the 
relativised quark model calculation agrees with the lattice results. 

The electric form factors calculated at q ~ 0.4 GeV for neutral decuplet baryons are 



illustrated in figure [13]. Neutral octet baryons are also given for reference. For the hyperons 
the form factors are dominated by the net charge of the light quarks. Octet and decuplet 
hyperons display similar behavior. For A the symmetric spin state of the quarks causes the 
form factor to vanish. This contrasts the neutron where the octet-spin asymmetry of the 
three quarks gives rise to a negative squared charge radius. 

With knowledge of quark distribution radii, the charge radii of neutral baryons may be 
calculated as in r 2 = J2i= s ,s,u e « r f f° r — °- Table [VI] summarizes the results. The lattice 
result for the squared charge radius of the neutron r^/r 2 = —0.11^°'^° encompasses the 



13 



experimental value |4l1|j4~T|| of —0.167(7). 



D. Magnetic-dipole form factors 

Our calculation of magnetic-dipole form factors is done at the smallest finite value of q 2 
available on our lattice. Table |VII| summarizes the form factors in units of natural magnetons 
(e/2Mg) where the mass of the baryon appears in the definition of the magneton. The 
magnetic moment \i is defined at q 2 = as £mi(0) = fi/(e/2M B ) and therefore we must 
scale our results from {?mi(<? 2 ) to ^mi(O). Lattice extrapolations in q 2 to q 2 = suffer 
from large statistical errors. To make contact with the experimental magnetic moments, we 
assume equal scaling of electric and magnetic form factors as in the octet baryon analysis 
and as supported by experimental measurements of nucleon form factors. In hyperons, the 
strange and light quark sectors are scaled separately by 

GmM g|(g) rs ^ 

G s mi{<?) GhWY 1 ■ ' 

and similarly for the light quarks, such that the magnetic moment of a hyperon is given by 

0in(O)=0W°)+%i(°). ( 3 - 4 ) 

where I labels the light quarks. For A baryons it is not necessary to separate the u- and 
ci-quark sectors due to the S'f/(2)-isospin symmetry of the correlation functions. Table 
|V111| gives the magnetic-dipole form factors for quarks of unit charge in units of natural 
magnetons. Combined with the results of table |T| the magnetic moments may be calculated. 

In the octet baryon analysis it was found that the magnitudes of the lattice results 
for magnetic moments were consistently smaller than the experimental measurements. It 
was argued that at f3 = 5.9 some deviations from asymptotic scaling may occur. A recent 
analysis |10| determines nucleon form factors at (3 — 6.0 on a cubic lattice with physical 



spatial dimensions roughly equal to our smaller y and z dimensions. Some improvement is 
seen in the magnitudes of the magnetic moments which are still 10 to 15% low compared 
to experiment. Chiral dressings of the nucleon may cause our linear extrapolation of the 
magnetic moments in l/n to underestimate the magnetic moments in the physical regime 



35] . Finite volume effects may also give rise to the underestimation of the magnetic moments 
as the baryon is restricted by its periodic images. The proton rms electric charge radius at 
% indicates the proton largely fills the lattice in our smaller y and z spatial dimensions. 
Non-quenched corrections may also provide additional contributions |^2| . 



To reduce the effects of these uncertainties, ratios of the lattice results to the lattice 
proton result are used when making comparisons with experimental measurements or model 



calculations. Table [IX| reports the extrapolated decuplet baryon magnetic moments in units 
of nuclear magnetons (/iAr), ratios of decuplet baryon moments and the proton magnetic 
moment and ratios scaled to reproduce the proton magnetic moment. These scaled ratios are 



illustrated in figure [14]. The expected qualitative behavior of mass dependence is displayed 
here. For example as d quarks are replaced by s quarks in going from A~ to S*~ through 
to Or the magnetic moments of the negatively charged hyperons decrease in magnitude. As 
in the case of the electric charge form factors of neutral baryons, the magnetic moments are 
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dominated by the net charges of the light quarks. S , f/(2)-isospin symmetry causes the A 
moment to vanish. 

The simple quark model formula for the magnetic moment of a decuplet baryon is simply 
given by the sum of the intrinsic moments of the quarks composing the baryon. Results using 
intrinsic moments determined by the p, n, and A moments |33j are also indicated in figure 



by horizontal dashed lines. The agreement is striking and suggests a baryon independence 
of the lattice effective quark moments. 

Figure [15] displays the lattice effective quark moments for quarks of unit charge within 
decuplet baryons. Effective moments have been defined by equating the lattice quark sec- 
tor contributions to the same sector of the S{7(6)-magnetic-moment formula derived from 
S'[/(6)-spin-flavor symmetry wave-functions. For example, S'f/(6)-spin-flavor symmetry gives 
the simple quark model formula 

fi p = A/j, u - Bfi d , (3.5) 

for the magnetic moment of the proton, where A = 4/3 and B = 1/3. The effective 
moment of the lattice w-quark in the proton is defined by equating the lattice w-quark sector 
contribution to |/x u , the corresponding w-quark sector contribution of the simple quark 
model. Similarly, the lattice cZ-quark sector contribution is equated with —\y>d m defining 
the effective (i-quark moment. In A + the lattice u-quark sector is equated with 2/j, u and the 
d-quark sector equals the effective <i-quark moment. Numerical values for effective quark 
moments are summarized in table |X| for both octet and decuplet baryons. 

Approximate baryon independence of decuplet baryon effective quark moments is dis- 
played in figure |l|. Closer examination of effective u-quark magnetic moment ratios does 
reveal some small baryon dependence of the quark moments. These results are summarized 
in table |X]]. The quenching of effective quark moments in hyperons is largely due to the 
baryon mass setting the scale at which quarks contribute to the baryon moment [|42|| . These 



results are in sharp contrast to the enormous baryon dependence of the effective quark 
moments seen in octet baryons as illustrated in figure |l^. 

It is instructive to consider ratios of the octet baryon effective quark moments with their 
decuplet partners. These ratios are illustrated in figure In the simple quark model all 
of these ratios are equal to 1. Dramatic effects are seen for the singly represented quarks 
u n , ss and v& where deviations from Si7(6)-octet wave-functions are large in the lattice 
results. All quarks show a significant baryon dependence with the exception of u p /u^ and 
possibly ug/us*. This rich octet /decuplet baryon dependence in the effective quark moments 
contrasts that of the quark charge distribution radii as illustrated in the ratios of figure |TT| . 

While lattice octet baryon moments differ significantly from the predictions of SU(6), it 
was not possible to determine how each of the coefficients A = 4/3 and B = 1/3 of ( ^.5| ) 
are altered. Using ratios of the u- and (i-quark sector contributions to p or n and isospin 
symmetry, the octet baryon analysis indicates 

| = 0.13«?<i (3.6) 

where 1/4 is the SU(Q) prediction. The decuplet baryon results indicate the effective 
moments and charge radii of doubly represented quarks are largely unchanged suggesting 
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A ~ 4/3 and indicating B is better approximated by 1/6 than the standard SU(6) value of 
1/3. 

Considerable effort has gone into model calculations of the A ++ and Q magnetic mo- 
ments. Here we collect together a handful of these calculations for comparison with our 
lattice results. Figure ^ displays ratios of A ++ jp magnetic moments. The model calcu- 
lations include results from the simple quark model |33| (Q.M.), cloudy bag jy| (C.B.), 
Skyrme (Skyr.), Bethe-Salpeter g§ (B.S.) and QCD sum rule g§ (S.R.) analyses. The 
experimental result (Expt.) is from a recent pion bremsstrahlung analysis |47j and there- 
fore has some model dependence. Earlier papers reported larger values [5^,55] in better 
agreement with our lattice result. 

Figure [19] displays ratios of Q~ jp magnetic moments. The labels are as in figure [18| with 
additions of an additive quark model [ED] in which effective masses (E.M.) of the quarks 
are used to estimate the intrinsic quark moments and a calculation in which relativistic 
corrections (R.C.) to baryon magnetic moments are considered [51]. Experimental results 
(Expt.) are from a recent investigation where Q~ hyperons are produced by a polarized 



neutral beam spin transfer reaction 52 



In the simple quark model the ratio of Q~ /A magnetic moments is 3. The lattice results 
suggest that spin-dependent forces give rise to a larger strange quark moment in Q~ relative 
to that in A as indicated in figure 17. The lattice ratio of Q~ /A magnetic moments is 



3.6^'° suggesting an enhancement over the simple quark model ratio. 



-0.6 



To complete the discussion of decuplet baryon magnetic-dipole properties figure £0 
summarizes our calculation of magnetic radii normalized by the magnetic moment as in 
J < r 2 > /Gmi(0)- The magnetic radii follow a similar pattern to that of the charge radii. 



E. Electric-quadrupole form factors 

The results of the higher-order multipole moments must be regarded as preliminary due 
to the lack of symmetry in the spatial dimensions of our lattice. The elongated x dimension 
coupled with possible finite volume restrictions may induce deformations from spherical 
symmetry. On the other hand, it is useful to examine these quantities with a view toward 
determining the feasibility of extracting E2 moments in future calculations. 

The E2 form factor is particularly interesting since it provides a glimpse into the shape 
of the baryon ground state. A small E2 form factor would cast serious doubt on Skyrme 
models where the hedgehog Skyrmion has an inherently large quadrupole moment. Our 



focus here is on estimating an upper bound for the E2 form factor. Tables |XIII| and |XIV 
summarize the lattice results for E2 form factors. 

The negative sign of the central values for positively charged baryons is consistent with 
deformation one might expect from our elongated lattice. Equation ( |2.23b| ) isolates the spin 
matrix element ( |s' | Sq 2 - 1 | |s ) and therefore determines the asymmetry of the ground state 
wave-function having some overlap with the spherical harmonic l^oW where 9 is measured 
relative to the spin quantization axis z. A negative E2 moment corresponds to an oblate 
shape which may be due to some broadening of the wave function in the longer x direction 
where overlap with periodic images is minimized. 

In nonrelativistic models where angular momentum I and spin s are constants of the 
motion, it is useful to consider angular momentum selection rules. Consider the transition 
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( 3 =1' + S'\ j-y \j=l + S ), 



(3.7) 



where j 1 is the total angular momentum of the photon and determines the multipole char- 
acter of the transition. Total angular momentum selection rules require j + j 1 = j' and 
therefore 

\j-f\<h<j+f. (3.8) 

Combined with parity conservation, ( |3.8p indicates there are at most 4 multipole form factors 
available in spin-3/2 to spin-3/2 transitions. Electric form factors correspond to s — s'. This 
yields the orbital angular momentum selection rule 



\l ~ l'\ <h<l + I', 



(3.9) 



indicating the E2 form factor vanishes in nonrelativistic models unless some configuration 
mixing is included in the baryon ground state. 

The relationship between our definition of the E2 form factor given in ( p. 16a ) and that 



usually calculated in models may be easily established by writing (|2.16a|) using the spherical 
harmonics in the Breit frame with s = s' = j = 3/2. Noting that 



[q x q] 



(2) 




*2m(?), 



(3.10) 



equation (|2.16a|) may be written as 
q 3 



( 1/(0) 



2' 2 



2 Att 

G E0 (q 2 ) + -rJ-Y 20 (q)G E2 (q 2 ) 



f.fiy<wi-f.i>- 



Model approaches determine current matrix elements via the fourier transform 

Using the small \q\ spherical harmonic expansion 

oo i 
*** = 47r E (2l + m & E Y *m(r) Y lm (q) , 

1=0 \ I" m=—l 

and choosing q = qz, it is straight forward to show 

g E2 {0) = M 2 B { d 3 r^(r) (3z 2 - r 2 ) ^(r), 



(3.11) 



(3.12) 



(3.13) 



(3.14) 



where 3z 2 — r 2 is the standard operator used for quadrupole moments. The factor M\ 
indicates that Q E 2^) gives the E2 form factor in units of (e/M§). 

A chiral quark model calculation gives Ge2{$) = —0.30 fm 2 for A ++ which lies 2a 
outside our bound of —0.03(13) fm 2 . However it may be more appropriate to compare the 
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dimensionless ratio of the E2 form factor and the squared charge radius. The generic Skyrme 
model result is 153IJ53 



/ 2 \ = ~^ J 3 » 3 - 15 

(r 2 ) v 25 

where ( r 2 ) y is the isovector charge radius and J3 is the 3-component of isospin. The 
analogous lattice ratio is 



^-^-^AT- (3-16) 



)y V 2; ( r 2)A 



The isospin dependence of the lattice results is proportional to the baryon charge which is 
not the case for the Skyrme model. For A ++ the lattice ratio is —0.08(30) which encompasses 
the Skyrme model ratio of —0.24. Future high statistics lattice calculations on a cubic lattice 
should be able to confirm or reject the hedgehog Skyrmion description of baryons. 



F. Magnetic-octupole form factors 

The magnetic-octupole form factors are large for heavy quark masses and appear to 
decrease as the quark mass becomes lighter. In the physical regime the uncertainties are 
sufficiently large to make the M3 moments consistent with zero with the exceptions of H*~ 
and Q~ hyperons. Tables |XV| and |XVI| summarize the lattice results for decuplet baryons 



and their quark contributions respectively. 

Angular momentum selection rules for M3 moments require s + 1 = s' limiting nonvan- 
ishing transitions to 

\\l-l'\-l\<j~,<l + l' + l- (3.17) 

Once again M3 transitions require nonzero orbital angular momentum admixtures in the 
ground state wave-function. Connection to model calculations of M3 moments may be 
made from our definition of the M3 moment in (|2.16b|) in a manner analogous to that for 



the E2 moment. M3 moments in hedgehog models 5IJ are 1/N C suppressed relative to Ml 
moments and cannot been determined using conventional semi-classical methods. 

The dimensionless ratio of interest for M3 moments relates the M3 form factor to the 
Ml form factor and the squared charge radius. For Q~ the ratio is 

LjM:i 1.4(4), (3.18) 



Gmi ( r 2 ) 

indicating a relatively large M3 moment. This ratio is much larger than one might expect 
from our discussion of the E2 form factor and angular momentum selection rules. For Q~ 
the ratio is 

^fy = 0.02(3). (3.19) 

However, in fully relativistic calculations angular momentum and spin on their own are no 
longer good quantum numbers. Gluons may also carry angular momentum allowing radical 
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changes from standard S'?7(6)-spin-flavor symmetry as seen in the magnetic properties of 



octet baryons and evidenced in the EMC polarized muon-proton scattering experiment [55 



Furthermore, lattice QCD is a relativistic theory with particle creation and annihilation al- 
lowing some overlap with mesonic dressings of baryons even in the quenched approximation. 

This pattern of small electric effects and large magnetic effects is reminiscent of the 
electromagnetic properties of octet baryons. The lattice results show a spin dependence in 
the quark distributions that accounts for the negative squared charge radius of the neutron 
by slightly increasing the <i-quark distribution radius relative to the w-quark radius. However, 
spin dependence in the magnetic properties is huge. For example, figure |T6] shows that the 
effective magnetic moment of a u quark in the neutron is roughly half that of the (i-quark 
when normalized to unit charge. 



IV. OVERVIEW AND OUTLOOK 

We have presented a fully relativistic formalism for isolating and extracting the four 
electromagnetic multipole form factors of spin-3/2 systems in lattice field theory. Results 
of the first lattice QCD analysis of 5'[/(3)-flavor decuplet baryons have been systematically 
examined to reveal new aspects of the underlying nonperturbative quark-gluon dynamics. 

The EO and Ml correlation functions for decuplet baryons show a broad plateau region 
allowing a reliable extraction of the electromagnetic form factors. Statistical uncertainties 
are similar to that seen in our octet baryon analysis. 

The qualitative mass dependence of decuplet baryon charge radii follows the anticipated 
pattern produced by a more localized strange quark distribution. The electric form factors 
of neutral baryons are dominated by the net charge of the light quarks as in the octet baryon 
results. Model calculations of the Q~ charge radius vary widely and the lattice results favor 



a calculation based on a relativised quark model |37] over the other bag and Skyrme models 
considered. 

Closer examination of the baryon charge radii reveals a behavior in the quark charge 
distribution radii consistent with a spin dependent force having an inverse relationship with 
the quark mass. This behavior is what one expects from the hyperfine interaction term 
of the one-gluon-exchange potential. The spin dependent force counteracts center-of-mass 
shifts and suppresses baryon dependence of the quark charge distribution radii. 

The center-of-mass shifts that give rise to an enhancement of the u-quark magnetic 
moment in the octet S° relative to that in the neutron are offset by this repulsive force in 
decuplet baryons. Variation in the quark effective magnetic moments from baryon to baryon 
is minimal for decuplet baryons. The residual effect is largely due to the baryon mass setting 
the scale at which quarks contribute to the baryon magnetic moment. 

The symmetric role of quarks in decuplet baryons contrasts that in octet baryons. The 
role of a quark in an octet baryon correlation function is different depending on whether the 
quark is singly or doubly represented. This difference, which cannot arise in the decuplet 
case, causes large variations in the electromagnetic properties of quarks in octet baryons . 

The lattice predictions of decuplet baryon magnetic moments agree with the simplest of 
quark models to a remarkable extent. This is in sharp contrast to the octet baryon analysis 
and is due to an approximate baryon independence of effective quark magnetic moments 
in decuplet baryons. The lattice prediction of the A ++ /p magnetic moment ratio is large 
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compared to a recent experimental analysis but compares better with previous analyses. 
Our Q~ magnetic moment agrees with experiment [52] and prefers the simple quark model 
result and a calculation based on effective quark masses |5(J over the other four model 
calculations considered. 

The E2 and M3 form factors depend on subtle differences in the correlation function 
ratios and as a result have larger statistical uncertainties. The M3 form factors are finite 
at all the values of k considered. The E2 form factor is at the threshold of confirming 
or rejecting the hedgehog Skyrmion description of baryons. A future high statistics lattice 
calculation on a cubic lattice will provide considerable insight. 

For electromagnetic form factors to provide a reliable test of QCD, one must understand 
and eliminate systematic errors. While the results of chiral perturbation theory may be used 
to assess the magnitude of possible systematic errors in the lattice extrapolation procedure 



35| , there can be no substitute to actual calculations probing regimes of lighter quark mass. 

A calculation of electromagnetic properties in full QCD would also assist in understanding 
systematic errors. However to calculate in full QCD, the contributions of disconnected quark 
loops as illustrated in figure [21] must be understood first. These loop contributions may 
play a key role in removing the discrepancies between lattice and experimental violations 
of magnetic moment sum rules | j42[| . Closed quark loop contributions may be estimated 
even in the quenched approximation and present a major challenge for future lattice QCD 
calculations. Such a calculation will provide insight into important questions such as the 
strangeness content of the proton and the general role played by sea quarks in hadrons. 

The examination of decuplet baryon structure through the lattice field theory approach 
to QCD has given a new and more detailed understanding of nonperturbative interactions. 
The underlying theme of the results presented here involves a cancellation of spin- dependent 
forces and center-of-mass effects which results in structure that is consistent with nondy- 
namical models such as the simplest quark model based on SU (6)-spin-flavor symmetry 
broken only by the quark masses. 

A future investigation |14[ will address Nj — > A transitions where there have been 
significant experimental efforts to measure the Ml and El transition moments. This, of 
course, has been accompanied by a plethora of model calculations for these quantities. The 
lattice results will be instrumental in assessing the reliability of the model analyses. 

The odd-parity spin-3/2 baryon octet also offers an interesting forum for checking model 
predictions. There has not been a lattice investigation of these states. The iV(1520) is the 
lowest lying baryon with I(J P ) — |(| ). It is expected to be stable with our present lattice 
parameters and therefore is accessible to lattice calculations. The neighboring (1535) is 
a source of possible contamination in the j = 3/2 interpolating field and therefore angular 
momentum projection operators may be required p6| . Calculations of the electromagnetic 



properties of iV| (1520) may be completely different from model expectations if gluons 
carry a significant fraction of the angular momentum usually attributed to quark degrees 
of freedom alone. With the anticipated experimental program at CEBAF, a lattice QCD 
analysis of N-^ 7 — ► transition moments is also of interest. 

It is very encouraging that lattice QCD evaluations of hadronic electromagnetic form 
factors give, for example, a pattern of magnetic moments which is as good or better than 
hadronic models with adjustable parameters which in some cases have been tuned to fit 
magnetic moments. With further study to understand and eliminate systematic errors, we 
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believe the study of hadron electromagnetic form factors will prove ultimately to be one of 
the best quantitative testing grounds for nonperturbative QCD. 
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TABLES 





TABLE I. 


Baryon Mass in 


lattice units (M#a). 




Baryon 


«i = 0.152 


k 2 = 0.154 


k 3 = 0.156 


K cr = 0.159 8(2) 


N 


1.09(3) 


0.96(3) 


0.84(3) 


0.61(5) 


A 


1.13(3) 


1.02(4) 


0.90(5) 


r\ r—T\ / rr\ 

0.70(7) 


S* 


1.13(3) 


1.05(4) 


0.98(4) 


0.84(5) 


"* 


1.13(3) 


1.09(4) 


1.05(4) 


0.98(4) 


n 


1.13(3) 






1.13(3) 




TARLF, TT 

1 ii !_/ J_J J_J i_ _1 . 


ud.i yuii cicLtiiL 


f narffp fnrrn TPtnirwa 

Lllal e,C 1U1 111 Idl^UVJlo. 




Baryon 


Ki = 0.152 


k 2 = 0.154 


k 3 = 0.156 


K cr = 0.159 8(2) 


A++ 


1.723(17) 


1.681(27) 


1.636(55) 


1.562(80) 


A+ 


0.861(9) 


0.841(13) 


0.818(27) 


0.781(40) 


A 


0.000 


0.000 


0.000 


0.000 


A — 


— 0.861(9 J 


—0.841(13) 


— 0.818(2 1 ) 


— U. (81(40 J 


£*+ 


0.861(9) 


0.831(13) 


0.796(28) 


0.743(40) 


S*° 


0.000 


-0.008(1) 


-0.017(5) 


-0.032(8) 


s *- 


-0.861(9) 


-0.847(11) 


-0.831(19) 


-0.805(28) 




0.000 


-0.016(3) 


-0.036(8) 


-0.065(14) 




-0.861(9) 


-0.854(10) 


-0.846(12) 


-0.832(16) 


n~ 


-0.861(9) 






-0.861(9) 


TABLE III. Electric charge form factors for single quarks of unit charge. 


Quark 


m = 0.152 


k 2 = 0.154 


k 3 = 0.156 


K cr = 0.159 8(2) 


UA 


0.861(9) 


0.841(13) 


0.818(27) 


0.781(40) 




0.861(9) 


0.839(12) 


0.814(24) 


0.774(34) 




0.861(9) 


0.863(9) 


0.866(13) 


0.869(18) 


■us* 


0.861(9) 


0.838(12) 


0.810(19) 


0.767(27) 




0.861(9) 


0.862(9) 


0.863(9) 


0.865(11) 


SQ 


0.861(9) 






0.861(9) 
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TABLE IV. Baryon rms charge radii normalized to unit charge in lattice units ( r 2 /a 2 ) 



Baryon 


ki = 0.152 


k 2 = 0.154 


K3 = 0.156 K cr 


= 0.1598(2) 


A++ 


3.71(13) 


4.02(19) 


4.34(39) 


4.90(57) 


A+ 


3.71(13) 


4.02(19) 


4.34(39) 


4.90(57) 


A~ 


3.71(13) 


4.02(19) 


4.34(39) 


4.90(57) 


£*+ 


3.71(13) 


4.15(19) 


4.63(38) 


5.42(56) 




3.71(13) 


3.92(16) 


4.16(28) 


4.54(41) 




o. / 11 ±o J 






/I 1 f^O^ 


o 


o 7i/iq\ 






o 71 / 1 
o. i l^±o ) 


TABLE V. 


Rms charge radii for 


single quarks 


of unit charge in lattice units ( r 2 


n . 1 /O 

/a 2 ) 1/2 . 


Quark 


«i = 0.152 


k 2 = 0.154 


k 3 = 0.156 K cr 


= 0.1598(2) 


UA 


3.71(13) 


4.02(19) 


4.34(39) 


4.90(57) 




3.71(13) 


4.03(18) 


4.40(33) 


4.99(48) 


SE* 


3.71(13) 


3.68(14) 


3.65(20) 


3.60(28) 


UE* 


3.71(13) 


4.05(17) 


A A A i c\ rr\ 

4.44(27) 


5.07(38) 


5"* 


3.71(13) 


3 70(13) 


3 68(14) 


3 66(16) 




3.71(13) 






3.71(13) 


TABLE VI 


. Neutral baryon rms 


charge radii 


in units of the proton radius ( |r 2 / 


n, .1/2 


Baryon 








Radius 


n 








0.34(22) 


A 








0.39(8) 


S° 








0.44(4) 


s° 








0.53(10) 










0.39(4) 










0.57(9) 
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TABLE VII. Baryon magnetic-dipole (Ml) form factors in units of natural magnetons 
{HB = e/2M B ). 



Baryon 


ll! = 0.152 


k 2 = 0.154 


k 3 = 0.156 


K cr = 0.159 8(2) 


A++ 


4.48(30) 


4.48(42) 


4.45(60) 


4.44(87) 


A+ 


2.24(15) 


2.24(21) 


2.22(30) 


2.22(44) 


A 


0.00 


0.00 


0.00 


0.00 


A — 

A 


—2.24(15) 


—2.24(21) 


— 2.22(30) 


— 2.22(44) 


£*+ 


2.24(15) 


2.36(21) 


2.47(30) 


2.69(44) 


S*° 


0.00 


0.058(13) 


0.119(45) 


0.224(75) 


Y,*- 


-2.24(15) 


-2.25(19) 


-2.24(24) 


-2.24(32) 




0.00 


0.119(23) 


0.235(68) 


0.45(12) 




-2.24(15) 


-2.25(17) 


-2.24(19) 


-2.24(23) 


n~ 


-2.24(15) 






-2.24(15) 


TABLE VIII. Magnetic-dipole (Ml) form factors for single quarks of unit charge in units of 


natural magnet 


ons (fi B = e/2M B ). 








Quark 


ki = 0.152 


k 2 = 0.154 


k s = 0.156 


K cr = 0.159 8(2) 


UA 


2.24(15) 


2.24(21) 


2.22(30) 


2.22(44) 




2.24(15) 


2.30(20) 


2.36(27) 


2.46(38) 


ss* 


2.24(15) 


2.13(17) 


2.00(20) 


1.79(25) 


"US* 


2.24(15) 


2.36(19) 


2.48(24) 


2.69(32) 




2.24(15) 


2.19(16) 


2.12(17) 


2.02(19) 


sn 


2.24(15) 






2.24(15) 


TABLE IX. 


Lattice results for decuplet baryon maj 


metic moments in units of nuclear magne- 


tons (//jv)j ratios of decuplet baryon 


moments to the proton magnetic moment and ratios scaled 


to reproduce the proton magnetic moment. 






Baryon 


Magnetic moment 


Ratio {n/ Hp) 


Scaled 










(mjv) 


A++ 


4.91(61) 




2.18(32) 


6.09(88) 


A+ 


2.46(31) 




1.09(16) 


3.05(44) 


A 


0.00 




0.00 


0.00 


A~ 


-2.46(31) 




-1.09(16) 


-3.05(44) 


£*+ 


2.55(26) 




1.13(14) 


3.16(40) 


S*° 


0.27(5) 




0.118(24) 


0.329(67) 


s *- 


-2.02(18) 




-0.90(10) 


-2.50(29) 




0.46(7) 




0.206(35) 


0.58(10) 




-1.68(12) 




-0.744(85) 


-2.08(24) 


Si- 


-1.40(10) 




-0.621(78) 


-1.73(22) 
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TABLE X. Effective quark magnetic moments for quarks of unit charge in units of nuclear 



magnetons (hn)- 





Octet Baryons 


Decuplet Baryons 


Quark 


Effective moment 


Quark 


Effective moment 


Up 


2.41(2 I ) 


ua 


2.46(31) 


dp 


1.10(74) 


d,A 


2.46(31) 




2.12(23) 


uz* 


2.29(22) 




0.59(37) 




1.49(13) 




1 53(271 


Us* 


2.14(18) 




1.24(10) 


s s * 


1.45(10) 


S A 


1.25(15) 


sn 


1.40(10) 




TART.F! XT Flffpptivp ?/- nnrl c-nnarlf 


magnetic moment ratios. 


Quark Ratio 






Lattice Result 


UY,*/UA 






0.93(5) 








0.94(3) 


us* /ua 






0.87(7) 








0.97(5) 


sn/s-,* 






0.97(3) 








0.94(7) 




TABLE XII. Baryon Magnetic Radii J ( , 


" 2 ) 1 {Quia 2 ) in 


lattice units. 


Bar yon 


«i = 0.152 k 2 = 0.154 


k 3 = 0.156 


K cr = 0.159 8(2) 


A++ 


71 H ^ 4 09(1 Q"l 


4.34(39) 


4.90(57) 


A+ 


3.71(13) 4.02(19) 


4.34(39) 


4.90(57) 


A~ 


3.71(13) 4.02(19) 


4.34(39) 


4.90(57) 


£*+ 


3.71(13) 4.13(19) 


4.58(38) 


5.34(56) 




3.71(13) 3.93(17) 


4.19(29) 


4.59(41) 




3.71(13) 3.82(14) 


3.98(18) 


4.22(23) 


n- 


3.71(13) 




3.71(13) 
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TABLE XIII. Baryon electric-quadrupole (E2) form factors in units of (e/M|). 



Baryon 


ki = 0.152 


k 2 = 0.154 


k 3 = 0.156 


K cr = 0.159 8(2) 


A++ 


-0.6(8) 


-0.8(14) 


-0.7(28) 


-1.0(39) 


A+ 


-0.3(4) 


-0.4(7) 


-0.4(14) 


-0.5(19) 


A 


0.0 


0.0 


0.0 


0.0 


A- 


0.3(4) 


0.4(7) 


0.4(14) 


0.5(19) 


£*+ 


-0.3(4) 


-0.3(7) 


-0.1(13) 


-0.1(19) 


s*° 


0.00 


0.05(7) 


0.15(23) 


0.26(36) 




0.3(4) 


0.4(6) 


0.4(9) 


0.5(13) 




0.00 


0.07(12) 


0.15(31) 


0.29(54) 




0.3(4) 


0.3(5) 


0.4(6) 


0.4(8) 


n- 


0.3(4) 






0.3(4) 


TABLE XIV. 


Electric-quadrupole (-E2) form factors for 


single quarks of unit charge in units 


of (e/M|). 










Quark 


«i = 0.152 


k 2 = 0.154 


k 3 = 0.156 


K cr = 0.159 8(2) 


«A 


-0.3(4) 


-0.4(7) 


-0.4(14) 






-0.3(4) 


-0.3(6) 


-0.3(11) 


-0.3(16) 




-0.3(4) 


-0.4(5) 


-0.7(6) 


-1.0(8) 




-0.3(4) 


-0.3(6) 


-0.2(9) 


-0.2(12) 




-0.3(4) 


-0.4(4) 


-0.5(5) 


-0.6(6) 


sn 


-0.3(4) 






-0.3(4) 


TABLE XV. Baryon maj 


metic-octupole (M3) form factors in 


units of (e/2M|). 


Baryon 


«i = 0.152 


k 2 = 0.154 


k 3 = 0.156 


K cr = 0.159 8(2) 


A++ 


113±30 


99±51 


80±73 


55±112 


A+ 


56±15 


50±25 


40±37 


27±56 


A 


0±0 


0±0 


0±0 


0±0 


A" 


-56±15 


-50±25 


-40±37 


-27±56 


£*+ 


56±15 


54±27 


48±48 


43±72 


s*° 


0±0 


0.8±2.9 


-0.4±7.7 


0.3±1.4 




-56±15 


-52±22 


-49±34 


-41±49 




0±0 


2±5 


2±14 


6±24 




-56±15 


-54±19 


-54±25 


-51±32 




-56±15 


-56±15 


-56±15 


-56±15 
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TABLE XVI. Magnetic-octupole (MS) form factors for single quarks of unit charge in units of 
(e/2M|). 

Quark n x = 0.152 k 2 = 0.154 k 3 = 0.156 n cr = 0.159 8(2) 

56±15 50±25 40±37 27±56 

u s * 56±15 53±24 48±41 42±61 

s E * 56±15 51±17 49±22 41±28 

u E * 56±15 57±23 56±37 57±55 

s s * 56±15 53±16 52±19 48±22 

s n 56±15 56±15 56±15 56±15 
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FIGURES 

FIG. 1. Octet and decuplet baryon masses. The nucleoli mass has been used to set the scale. 
Horizontal lines indicate the experimental values. 

FIG. 2. Quark contributions to the correlation function ratios isolating the S* + electric-charge 
form factor as in ( |2.23a| ). The correlation functions are for k u = K3. Ratios for larger quark mass 
have smaller statistical uncertainties. 

FIG. 3. Quark contributions to the correlation function ratios isolating the H*° mag- 



netic-dipole form factor as in (2.23c). For purposes of illustration the leading kinematical factor, 
which varies in the second order jackknife procedure, has been set to a constant. The correlation 
functions are for k u = K3. Ratios for larger quark mass have smaller statistical uncertainties. 



FIG. 4. A comparison of the ratio sums in the parentheses of equations ( p. 23a ) and ( p. 23b ) 
isolating Geo and Qe2 respectively. Correlation function ratios are for k u = m. 



FIG. 5. A comparison of the ratio sums in the parentheses of equations ( |2.23cj ) and ( p. 23d ) 



isolating Qmi and Qm3 respectively. The M3 form factor is finite. However, the ratio fails to form 
a clear plateau. Correlation function ratios are for k u = k\. 

FIG. 6. Electric charge radii for charged octet and decuplet baryons in units of the proton 
charge radius. Radii of A ++ and A~ are equal to that of A + when normalized by the total baryon 
charge. 

FIG. 7. Extrapolations of the electric charge distributions of quarks within the octet S + 
baryon. The radius of the s-quark distribution decreases as the u quarks become lighter. 

FIG. 8. Extrapolations of the electric charge distributions of quarks within the decuplet S* + 
baryon. The radius of the s-quark distribution is largely unaffected by changes in the u-quark 
mass. 

FIG. 9. Charge distribution radii of unit charge quarks within decuplet baryons. Symmetric 
isospin symmetry equates u- and ci-quark properties in decuplet baryons. No significant baryon 
dependence of the charge radii is seen. 

FIG. 10. Charge distribution radii of unit charge quarks within octet baryons. Center-of-mass 
shifts and spin-dependent forces give rise to significant baryon dependence of the charge radii. 

FIG. 11. Ratios of octet and decuplet quark charge distribution radii. Doubly represented 
quarks show little octet/decuplet baryon dependence. 
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FIG. 12. The f2~ /p charge radius ratio. The model results include a calculation based on a 
relativised quark model (Q.M.), MIT bag (MIT) and cloudy bag (C.B.) models, and a bound state 
approach Skyrme (Skyr.) model. References may be found in the text. 

FIG. 13. Electric charge form factors of neutral baryons calculated at q ~ 0.4 GeV. The light 
quark charges dominate the hyperon form factors while decuplet spin symmetry causes the A 
form factor to vanish. 

FIG. 14. Ratios of the lattice baryon magnetic moments to the lattice proton result scaled to 
reproduce the proton magnetic moment. Horizontal dashed lines indicate the simple quark model 
predictions. 

FIG. 15. Effective quark moments in decuplet baryons. Approximate baryon independence of 
the quark moments is displayed. 

FIG. 16. Effective quark moments in octet baryons. Significant baryon dependence of the 
quark moments is illustrated. 

FIG. 17. Ratios of octet baryon effective quark moments with their decuplet partners. 
Octet/decuplet dependence is seen throughout the quarks with the exception of u p /ua and possibly 

UJ]/UJ]*. 

FIG. 18. Comparison of the lattice A ++ /p magnetic moment ratio (Latt.) with model calcula- 
tions and an experimental result (Expt.) having some model dependence. The model calculations 
are the simple quark model (Q.M.), cloudy bag (C.B.), Skyrme (Skyr.), Bethe-Salpeter (B.S.) and 
QCD sum rules (S.R.). References are given in the text. 

FIG. 19. Comparison of the lattice Q~ /p magnetic moment ratio (Latt.) with model calcu- 
lations and experiment (Expt.). The model calculations are as in figure |l8| with additions of an 
additive quark model based on effective quark masses (E.M.) and a calculation including relativistic 
corrections (R.C.). References are given in the text. 

FIG. 20. Magnetic radii of decuplet baryons in units of the proton magnetic radius. The 
magnetic radii follow a similar pattern to that of the charge radii. 

FIG. 21. A skeleton diagram of a disconnected quark loop contributing to the magnetic moment 
of a baryon. The diagram may be dressed with an arbitrary number of gluons. 



31 



